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Q . Abstract 

We present a numerical study of an SU(3) gauged 2D model for ad- 
joint scalar fields, defined by dimensional reduction of pure gauge QCD 
in (2+l)D at high temperature. We show that the correlations between 
Polyakov loops are saturated by two colourless bound states, respectively 
r-| ! even and odd under the Z<i symmetry related to time reversal in the orig- 

inal theory. Their contributions (poles) in correlation functions of local 
composite operators A n respectively of degree n = 2p and 2p + 1 in the 
5_i ■ scalar fields (p = 1, 2) fulfill factorization. The contributions of two par- 

ticle states (cuts) are detected. Their size agrees with estimates based on 
a meanfield-like decomposition of the p = 2 operators into polynomials 
in p — 1 operators. In contrast to the naive picture of Debye screening, 
no sizable signal in any A n correlation can be attributed to 1/n times 
a Debye screening length associated with n elementary fields. These re- 
sults are quantitatively consistent with the picture of scalar "matter" 
fields confined within colourless boundstates whose residual "strong" in- 
teractions are very weak. 
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1 Introduction 



Dimensional reduction is a powerful technique to study the infrared region of field 
theories at high temperature [0-0. Combined with a non-perturbative lattice sim- 
ulation of the reduced model, it has been employed to investigate the properties of 
gauge theories and QCD with dynamical quarks in the plasma phase 0-0]. For a 
recent review see [[Hj . 

It is, however, still not clear, what are the limitations and the domain of va- 
lidity of this approach. Therefore, in a recent work []TTJ, we have studied in detail 
the reduction to two dimensions of pure gauge QCD in (2+1) dimensions at high 
temperature. We refer the reader to this article for a more complete discussion of 
our motivations, references to the related literature and details on the reduction 
procedure. The reduced model is a model for scalars belonging to the SU(3) algebra 
(formerly the electric gluons in a static gauge). They interact with the 2D gauge 
fields, the static parts of the original 3D spatial gauge fields, and via an effective 
potential whose self-couplings are computed by perturbative integration over the 
non-static degrees of freedom. In [|ITJ we restricted the perturbative integration to 
the one loop order. The main conclusion of our investigation is that dimensional re- 
duction works very well in this case. In particular it was shown that it works within 
a few percent for the correlation function of Polyakov loops (as well as for spacelike 
Wilson loops) down to 1.5T C , where T c is the critical temperature in the (2 + 1)D 
theory. In fact for T > 1.5T C numerical simulations showed that the Polyakov loops 
correlations measured in (2+l)D QCD and in the reduced model were very close to 
one another down to quite short distances, although our formalism is in principle an 
expansion both in high temperature and in p/T, where p is the relevant momentum 
scale. It was further shown, by simulations at different values of the bare parame- 
ters, that our measurements are in the scaling region, and thus our results are valid 
in the continuum limit. 

Less expected was the observation that these correlations assume a shape typical 
of single particle propagation, as opposed to the standard picture of a screening mass 
associated with two massive electric gluons. This feature, together with our previous 
motivations, invites us to pursue the numerical exploration of the reduced model 
per se, in particular the investigation of states connected with the Polyakov loop 
correlations. 

While in JTTJ] we only measured the correlations of Polyakov loops, here we 
analyze separately those of SU(3) invariant polynomials of degree n in the elementary 
scalars A, namely A n = tr A n . The action is invariant under a global sign reversal of 
all the A's. This ^-symmetry, which we denote by R T , following the authors of Ref. 
12| , corresponds to Euclidean time reversal in the (2 + 1)D theory. In this latter 



article it was suggested to investigate operators odd under this symmetry, to obtain a 
possible gauge invariant definition of the Debye mass. To investigate both operators 
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which are even and odd under R T we consider separately correlations involving even 
and odd polynomials, n = 2p and 2p+l. The R T -symmetry may be spontaneously 
broken in the 2D model. In fact there exists two phases, corresponding to R T being 
conserved or broken. Only in the symmetric phase, the model corresponds to the 



reduction of the high temperature (2 + 1) QCD phase |T3[ , [[LJ] , |T2[ , |J . The details 



of the actual phase diagram will not be studied in the present article. From our 
data we can conclude, however, that the values of the coupling constants in the 
reduced 2D theory are in the unphysical broken phase. In the same way as has 
been done for the reduction from (3 + 1)D to 3D, we solve this problem by working 
in the metastable part of the symmetric phase. Using zero field initial conditions 
on large enough lattices, we make sure that we stay in the phase of unbroken R T , 
where even and odd operators do not mix. Investigations of states in the full and 
reduced model in the case of the (3+1) — ^3 reduction can be found in 0-0 and 
I5|— [20|, although a detailed analysis of the nature of these states, as we perform 



in this paper, have not yet been made. 

The 2D action under consideration is recalled in Section 2, together with its 
meaning in terms of the (2+l)D QCD model, from which it originates, and the 
relevant operators and correlations are defined. The simulations are performed in 
the temperature range 2T C to 12T C , where T c is the deconfining temperature of the 
latter model. Our results are presented and discussed in Sections 3 and 4. In Section 
3, we first show that the measured correlations fulfill the factorization properties 
expected if the lowest states in the n— even and n— odd channels are two distinct 
one particle states, whose masses are then extracted from the large distance decays. 



According to the criteria proposed in Ref.[12], the state found in the odd channel 



is a candidate to define a Debye screening mass, although not the only one. In 
Section 4, we further analyze the composite operators A n , n = 2p and 2p + 1 and 
their correlations A niTn , showing that all the properties observed for p = 2 can be 
deduced with a good accuracy from their knowledge for p = 1. This follows from 
the assumption that, given the SU(3) and R T symmetry constraints, the effective 
model for the elementary A's after integration over the gauge fields is a free field 
model for the massive composites A 2 and A 3 . In particular, we give evidence that 
the (small) deviations from factorization observed at short distances are mainly due 
to intermediate states containing two of the above particles. The summary and the 
conclusions can be found in a last Section 5. 
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2 The Reduced Action, Operators and Correla- 
tions 



In this section, we write down the reduced 2D action derived in [11|, and define our 
notations and the quantities of interest for the present work. 

The lattice is an L s x L s square; the spacing is a, set to one unless specified 
otherwise. The weight in the partition function is written exp(— S), with S a function 
of the SU(3) group elements U(x; i), % = 1, 2 on the links and of the scalars A(x) in 
the adjoint representation on the sites: 

8 1 
A(x) = ]TA Q 0r)A Q , tvX a X l3 ^-S a0 . (1) 

a=l 

Greek superscripts on A will always denote colour indices, unlike integers n, m used 
in powers of the algebra element A. We write the 2D reduced action as follows: 

S = Su + Su,a + Sa, (2) 
S v = /3 3 L V(l - -3ttTU(x] l)U(x + al;2)U(x + a2; l) _1 ^(x; 2)" 1 " 

2 



U,A 



(3 3 L 



^^tr(A(f/)A(x)) , (3) 



6 

X 1=1 



Di(U)A(x) = U{x; i)A(x + ai)U(x; z)" 1 - A(x), 
S A = J2 k 2trA 2 {x) + k 4 (tiA 2 {x)) 2 . 

X 

In the above, Su is the pure gauge term, Su,a the gauge invariant kinetic term 
for the scalars and Sa the scalar potential, whose self couplings k 2 and fc 4 result from 
the one-loop integration over the non-static components of the 3D gauge fields. All 
terms have the global i? T -symmetry A(x) — ► —A(x), while the Z3 symmetry of the 
original (2 + 1)D ST/ (3) model is broken by the perturbative reduction procedure. 
It was found in IITTI that 



3 / \ 5 

— [c Q \ogL Q + cij; c = 1, ci = -log2-l, (4) 



64tt 

The values of the parameters fa, L follow from the original lattice regularization of 
3D pure QCD at temperature T and gauge coupling squared g\ in the continuum. 
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The latter has dimension one in energy and is used to set the scale: 



l = -L. 

ai 

Accordingly, the continuum limit a — > is obtained by letting /3 3 and L go to 
infinity with the dimensionless temperature 



being kept fixed. 

The original three dimensional pure SU (3) gauge theory has a global Z 3 symme- 
try. The corresponding order parameter is the Polyakov loop. It is a static operator. 
In the reduced theory it has the form 

L{x) = -tr exp[iLoA(x)]. (7) 

At sufficiently high temperature the symmetry is spontaneously broken, signalling 
the deconfinement of static charges in the fundamental representation. The reduced 
theory in the above form does not have the Z 3 symmetry any more, because the 
perturbative reduction is made around one of the broken vacua, where A(x) = 0. 
The phase transition in the three dimensional theory appears at r c ~ 0.61 p3|. 
The reduced model should be valid in the deconfined phase, at sufficiently high 
temperature and long distances. In |ll| we employed it to investigate the correlations 



between Polyakov loops in this phase, where they are related to screening. We 
performed a detailed numerical analysis in the reduced model and a comparison 
with the results in the full (2 + l)D theory. Our simulations were performed for two 
values of the parameter L , namely 4 and 8. It was shown that scaling was very 
good, when comparing the data collected for fixed r at these two values of L . In 
this paper, we will show data collected for L = 4. As this is in the scaling domain, 
we can give the results in physical units. Distances R are given by 



RT = r/L, 



o- 



where r is the distance in lattice units, and temperatures are measured in units of 
the three dimensional critical temperature T c . For Lq fixed in the scaling region one 
may use 

T/T c = p 3 /p 3c . (9) 
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To discuss the effective Lagrangian of the model at fixed r in the scaling limit, 
it is convenient to normalize the scalar fields differently, defining <p{x) via 



A{x) = <j>{x) 



L0P3 



(10) 



The corresponding Lagrangian £ e // was derived in JTT| from the small a expansion 
of the effective action S. For clarity of the discussion we reproduce it here below. 
In S, we define Ai by 

U(x;i) = exp[iag2Ai(x)], (11) 

where g\ = g\T is the effective coupling of the 2D theory. Taking the limit a — > 
(but for the UV logarithm logLo = — log aT in the <fi 2 term), one obtains 



C 



eff 



D, 



F 
1 y 



C 



CT 



c=l 

9 i (f> + ig2[A i ,(f>), 
diAj - djAi + ig 2 [Ai, Aj], 
3^1 



(12) 



2tt 



-log(oT) + 5/2 log 2- 1 



tr 



(13) 



This is a 2D, SU(3) gauge invariant Lagrangian for an adjoint scalar 0, but it 
is far from being the most general one. The gauge coupling g 2 , with its canonical 
dimension one in energy, sets the scale. The non kinetic quadratic term is the coun- 
terterm Cct, suited to a lattice UV regularization with spacing a. The appearance 
of this term in the context of dimensional reduction in e.g. the lattice regularization 
framework was first discussed in Refs.|Q|[^]. It is well known that such logarithmic 
terms also appear in general, when one wants to define the continuum limit of a 
2D lattice model, see e.g.[^TJ|^2|. The reflection symmetry — > —0 is related to 
the euclidean time reflection in the original (2 + 1)D theory, noted R T and dis- 
cussed by Arnold and Yaffe JT2" |. In the two dimensional model of Eq. (^), it may, 
however, be spontaneously broken in some subspace of the unrestricted parameter 
space /33L0, k 2 , k^. As was discussed in Refs. [JT3[] (|TJJ 0, and as follows from the 
invariance of the original three dimensional theory under euclidean time reflection, 
the physical phase is the R T symmetric phase. The phase diagram in the 3D adjoint 
Higgs model, related to 4 — > 3 QCD reduction has been studied in Ref. [IB]. For the 
case of SU(2) in (3+l)D, R T is the center of the gauge group so that R T breaking 
is also gauge symmetry breaking, a subject previously discussed in Refs. fI5|| , p4 
and @. 

We have made a numerical investigation of the relative positions of the phase 
transition and of the reduction point in the reduced model, for T/T c = 1.97, and 
for two values L = 4 and 8. We find that the reduction point is near to the phase 
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transition, but in the broken phase. The phase transition is first order and strong 
enough, so that we can study the reduction in the metastable symmetric phase, 
by using appropriate starting values for the fields, and employing a large enough 
(32 x 32) lattice. 

We now turn to the definitions of the quantities relevant for the present investi- 
gation of the 2D model. The Polyakov loop correlation V(x) is defined by 

V{x) = (L(0)lS{x))-\(L)\*, (14) 

where L(x) is the Polyakov loop operator of Eq. (|7]). We now expand L(x) in powers 
of A(x), and we will study the operators A n (x) and connected correlations of their 
traces, defining 

A n (x) = trA n (x), (15) 
A n>m (x) = (A n (x)A m (0)) - (An) (Am). (16) 

When not ambiguous, the notation A n may also represent (A n (x)). Any operator 
A n {x) is gauge invariant, even or odd under the i? T -symmetry of S for n even or 
odd respectively. Because the reduced model was derived from a small A-fields 
expansion, its properties are significant for the (2+l)D model in the unbroken R T 
phase only, where Azp+i = and Ai v is small. 

In this article we will concentrate on the study of correlations of these operators, 
which are directly related to the static operator L(x) and therefore relevant for the 
high temperaure properties of the original (2 + 1)D SU(3) theory. As can be easily 
seen from the effective Lagrangian above, the 2D model has further symmetries 
beside R T , namely reflections of the 1- and 2-axis, which can be used to classify 
further operators, whose correlations may be studied. (In 2D there is no spin quan- 



tum number.) A corresponding analysis has been made in [2(J in the case of 
the three dimensional adjoint Higgs model. Although a full numerical analysis of 
the two dimensional adjoint Higgs model certainly has an interest per se, the oper- 
ators other than those which we defined above are not directly related to a static 
(2 + 1)D operator. One would need a further study to acertain to what extent their 
correlations are related to the high temperature physics of the original theory. We 
therefore do not discuss them in the context of this paper. 

We have performed a numerical simulation of the model defined above, with a 
flat measure for the A's and the standard De Haar measure for the gauge fields. 



The algorithm and error estimate techniques used are the same as in |IT] and not 
reproduced here. The lattice size is Ls = 32, and L = 4, throughout the present 
work. The (3% values are 29, 42, 84 and 173. This corresponds to values of T/T c 
equal approximately to 1.97, 2.85, 5.70 and 11.73 respectively. We were able to 
extract information from operators and correlations corresponding to n = 2 to 5. 



The cases n = 2 and 3 for the 3D reduced model were investigated in |TB . 
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3 The Lowest States of the Scalar Spectrum 



In this section we present our results for the A nrn correlations measured, and de- 
scribe them for each temperature in terms of two states S and P, respectively even 
and odd under R T and appearing for for n, m both even and both odd. Their phys- 
ical masses will be denoted M$ and Mp. For simplicity we often use in this and the 
next section the bare parameters r and /3 3 , related to R and T by Eqs. (|8]) and (|9|). 

In the simulations reported here, all runs were initialized with zero A-fields values 
and, as already stated, the system stayed in the metastable R T unbroken phase, as 
desired, with (A 2p +i(x)) being always compatible with zero. In this respect, the 
situation is thus similar to that encountered in ( [p~3|] - [fl6[] ) . 



For a first look at the data obtained in even and odd channels, we show in Fig. 1 
the on-axis correlations A njTn (r), n < m e [2, 5] at T/T c = 1.97 (/3 3 = 29). They are 
plotted against RT, that is the physical distance in units of the inverse temperature. 
In the even cases, the three correlations all have the same shape, and they decrease 
by about one order of magnitude each time two more powers of A are involved. The 
same is true for the odd cases, with a common decay of the correlations steeper than 
in the former case (smaller correlation length in lattice units). The overall situation 
is similar for /? 3 higher. For n, m larger than 5, as well as for /3 3 very large, the 
signal/noise ratio becomes very small. This can be understood at the qualitative 
level by noting that the rescaling Eq. ( |10|) of the A- fields normalizes the kinetic term 
for the 0-fields to the standard, parameter independent form 1/2 tr (Z}j0) 2 . Hence 
if the field renormalization by the interactions is weak, the <fi correlations should 
depend only weakly on /3 3 , which means that A n>m scales like p~\ n ^~ m >l 2 _ This will 
be illustrated more quantitatively in section 4. Due to this scale factor, the A- 
fields remain "small" in practice down to quite low values of /? 3 , which a posteriori 
explains why the perturbative reduction may still work at a temperature as low as 
1.5 T c . In fact, we checked that the Polyakov loop correlations are actually fully 
reconstructed within errors by keeping {n,m} up to {5,5} only in their expansion 
in A n m 's obtained from the small A expansion of (|T4"D. 

Now we want to analyze these A n ^ m data quantitatively in terms of the lowest 
states of the spectrum. Let rrii = a Mi be the lowest mass with quantum number 
i = S,P, in lattice units. For particle i we introduce a lattice propagator A Latt (m„ p) 
in momentum space: 



^lltti^P) = f + 4sinh(m 2 /4), (17) 
p 2 = 4sin 2 (p 1 /2) +4sin 2 (p 2 /2). 
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The corresponding contribution to A n ^ m (r) then reads 

A n , m (m u r) = shmjj cos (Pi r )^Latt{mi,p), (18) 

S Pl,P2 

where g l nm measures the residue of A n>m at the pole of QTTD, which on large enough 
lattices sits at p 2 ~ p 2 ~ — mf . With our definitions, is non zero only for i — S 
if n and m are even, and for i = P if n and m are odd. Prior to any fit of the masses 
to the data, we notice that a first consequence of our expectations on the lowest part 
of the spectrum comes from residue factorization g l nm = 7„7m, a property which we 
can probe directly on the correlations since, as r becomes large, it implies 



A ntn (r) A n+ 2 !n+2 (r) ^ 



A « = 72 (I) > l - l 19 J 



That it is so is demonstrated on Figs. (2-5) showing X 2 and X 3 (symbols o) for 
T/T c = 1.97 and 5.7 ((3^ = 29 and 84). The agreement is very good in all cases, 
although the quality of the data is poorer for X 3 due to the correlations involving 
A 5 getting very small. Similar results are obtained for other values of T/T c . We 
thus conclude at this point that single particle propagation accounts very well for 
the largest correlation length occuring in each of the two channels. Most of the 
observed deviations of X n from one will be interpreted in the next section in terms 
of two particle state contributions (symbols o in the same figures). 

We now proceed to assign values to the two lowest masses Ms and Mp expected 
from the above findings. This we do by various ways, in order to further enforce the 
statement that the correlations do have the characteristics associated with the pole 
structure of Eq. ([I7|). Down to r ~ 1, an excellent approximation to the on-axis 
correlation ([TJD is given by 



A n>m {m i ,r)~c( T \-e- m ^ + - 1 — e'^'A, (20) 
\[mir\y 2 [mi(L s - r)\ 1 ' 2 J 

where c is constant in r. We performed fits of this formula to all our A Ujm (r) data 
taken at r > r min . These fits are stable with respect to r min provided it is larger than 
about 4, and the values found for m ; in different correlations are always consistent 
with each other. The smallest errors were obtained by using fits to A 2t 2 and A 3 3 . 

Effective masses can also be obtained without any fitting by using 0— momentum 
correlations, defined for a generic x— space correlation C(xx,X2) by 

C°(r) = ^J2C(r,x 2 ). (21) 

X2 
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If the lowest mass in C is m, the ansatz ( |TTD gives 



C°(r) oc cosh(m(L s /2-r)Y (22) 

and m can be extracted at any r by inverting this relation: 

log(V(r) + V^ 2 W-l), (23) 
C°(r + 1) + C°(r - 1) 



Y(r) 



2C°(r) 



As an overall consistency check, we have extracted an effective mass m e ^(r) from 
0— momentum Polyakov loops correlations (|I3D, and compared it to the ms values 
obtained by our fits to A 2 2 . We find that m e ^ (r) is indeed nearly constant, in fact 
slowly decreasing towards a value compatible with ms, due to smaller and steeper 
contributions to fll4|) of the heavier particle P. 

A contrario, we invalidate the interpretation of the largest correlation length 
in v4 nn as being n times shorter than the "Debye screening length", the inverse 
of a mass m^ associated with "electric" gluons of the initial (2+l)D model (the 
scalars of the reduced model). This scenario was advocated by D'Hoker in his 
perturbative study of QCD 3 at high temperature [|4j]. If such was the case, the 
on-axis correlations A n>n should rather look like 

A „ JnmE , r) K (j^pe— + [mE(L \ r)]1/ ^ -- r >y, (24) 



which differs in shape from (EH), as was illustrated in O] for Polyakov loop cor 



relations. We nevertheless tried fits with fl24|) , but got definitely worse agreement, 
in the range of temperatures, which we have investigated, i.e up to 12T C . Hence 
the above picture is ruled out by the data in this temperature range, and if a mass 
can be defined for the electric gluon in high temperature QCD3 it is most probably 
larger than both m.5/2 and mp/3. In a "constituent gluon" picture, as advocated 
in Ref. [EJ, one would have bound states instead of a cut. One would, however, 
expect mp/ms ~ 3/2. 

Our final results for the S and P masses in units of the scale y ' g\T are collected 
in Table [I] for the values of T/T c investigated. They are taken from fits to A 2t 2 
and A3 3 respectively. The values for Ms agree with those obtained in [|ll]] from the 
Polyakov loop correlations. As can be seen from the tables, the ratios Mp/Ms vary 
with T/T c . There is, however, no clear tendancy in the region we have investigated, 
the ratios being 1.8, 2.0, 1.7, 1.6 in order of increasing temperature. We can, of 
course, not exclude that the ratio goes to 1.5 at still higher temperatures. 
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T/T c 


Ms/^MT 


M P /^fT 


1.97 


1.39(3) 


2.51(13) 


2.85 


1.48(5) 


3.02(16) 


5.70 


1.89(5) 


3.22(15) 


11.73 


2.12(8) 


3.44(11) 



Table 1: Masses in units of y g%T for the 5* and P states, as measured from fits to 
A 22 and A 3 3 respectively, for different values of T/T c . 

4 Weak "Strong" Interactions Between Colour- 
less States 

Here we will show that even at quite short distances (r small compared to m^ 1 ) 
all the condensates A n = (A n (x)) and correlations A n)m {r) can be reconstructed to 
a good accuracy from the data for A 2 , A 2)2 and A$$. The assumption is that the 
elementary fields A a (x) (Greek superscripts are colour indices) interact only through 
5* and P exchanges between the non-interacting composite A 2 (x) and A 3 (x), the 
scale of the fields being fixed by the size of A 2 , while A 3 =0. The precise way how 
this idea is implemented and the corresponding technicalities are detailed in the 
appendix. 

Here we limit ourselves to the simplest applications and give the results, starting 
with the local condensates. 



4.1 Weak Residual Interactions: The ^-fields condensates 

Since SU(3) has rank 2, any A n (x) can be reduced to a polynomial in A 2 (x) and 
A 3 (x). An elegant method p6 and explicit formulae are given in the Appendix. 
For n odd A n is zero by R T symmetry. For n even, we apply Wick contraction to 
all pairs of A a elementary fields, followed by the meanfield-like substitution 

A a (x)A?(x) -> \5 a , p A 2 (x). (25) 
As an illustration, consider A4. With the definitions of Section 2, we have 

1 8 

2A 4 (x) = A\{x) = 7^ A a {x)A a {x)A p {x)A p {x). (26) 

a,/3=l 

There we apply (p5|) and then replace A 2 (x) by its average A 2 . The A a A a and 
A 13 A 13 contractions give (8 x A 2 /4) 2 , and the additional contributions from a = (3 
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give 2 x 8(A 2 /4) 2 . Noting that A 2 , 2 (0) = (A%(x))-A$ (see (0)), the net result can 
be put into the two equivalent forms 

2A, = ^-Al (27) 

A 2 M = \a\. (28) 

This prediction is remarkably well verified in all cases. At /3 3 = 29, the left and 
right hand sides of (^8]) are respectively 4.86(1) 10 -3 and 4.825(10) 10~ 3 . They are 
6.093(13) 10~ 4 and 6.069(1) 10~ 4 at (3% = 84. Similar manipulations lead to 

A 3 , 3 (0) = ^Al (29) 

In this case the left and right hand sides are measured to be 2.105(7) 10~ 4 and 
2.095(6) 10~ 4 for /3 3 = 29, 9.365(30) 10~ 6 and 9.344(3) 10" 6 for (3 3 = 84. Hence the 
effects of residual interactions via non quadratic effective couplings in A 2 {x) and 
As(x) are less than the percent in the correlations at zero distance. 

Before going to the correlations at non zero distance, let us discuss their normal- 
ization, as measured by the values of A 2i2 (0) and A 3j3 (0) described just above. At 
the beginning of section 3, we argued that the behaviour in /3 3 , n, m observed for 
A n , m could follow from the absence of a large renormalization, by the interactions, 
of the 0-fields defined by Eq. ([TDD . Here we note that in the confined phase the 



effective degrees of freedom are the massive composites fa = trfa, i = 2, 3, so that 
in the limit where they are considered as free fields, one may write (see (|i~7|) ) 



(fa(O)fa(O)) ~ ft [d 2 p ^ \ (30) 
J p z + 4 smh(m i z /4) 

p 2 = 4sin 2 G» 1 /2) +4sin 2 Go 2 /2), 

the residue Ri being one if neither nor the composites get renormalized. We 
computed Ri as the ratio of the l.h.s. of (p0|), directly measured, to the integral 
in the r.h.s, evaluated numerically on the lattice for the mass values fitted to the 
correlation data. The result is shown in Fig.6: in the whole temperature range, both 
residues in the even and odd channels remain uniformly very close to one. 



4.2 Weak Residual Interactions: Properties of the Correla- 
tions 

As we have seen in section 2 (see Fig. 1), the different A niTO (r)'s corresponding to 
the same channel have very similar shapes. Their analysis in terms of one particle 
exchange was successful, confirmed by the agreement with residue factorization. 
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Nevertheless although the quantity X n of Eq. ( |i~9|) does go to one at large distances, 
it is significantly different from one at medium and short distances (see Figs. 2-5). 
We will now show that two particle exchange is responsible for most of this lack of 
factorization. 

The simplest consequence of our assumptions for correlations at finite r is, using 

Eq. 

A a>4 (r) = jA 2 A 2i2 (r), (31) 

which is very well verified at any distance as shown in Fig.7 for T/T c = 1.97 
(/?3 = 29). A similar agreement is found for the relation 

A 3 , 5 (r) = -A 2 A 3 , 3 (r), (32) 

derived in the appendix. A new situation arises when we consider A^ or ^5 
where both the initial and final states may couple to a two particle state, (S S) or 
(S P) respectively. Then the intermediate state in a connected correlation between 
and r may consist of either one or two particles. For example, to compute A^(r), 
we apply the substitution rule ( p5|) to the sum (|26|), and then average using the 
definitions of A 2 2 and A 2 . One finds 



5\ 2 



4AiA 2i2 (r) +2Ai 2 ( 



(33) 



A similar treatment given in the appendix leads to the prediction 

/3 c i\ 2 r 1 
A, 5 (r) = (-) [A 2 , A 3i3 (r) + A 2 , 2 (r) A 3 , 3 (r)\ . (34) 

In the two expressions above, the second contribution, a product of two propagators 
in space, is that of a two-particle intermediate state, and it provides a correction to 
exact factorization. From the definitions Eq.(|I^) of X 2 and X 3 , one actually gets 
the following estimates: 

X 2 (r) ~ X 2 (r) = l + ^g|i, (35) 
X 3 (r) ~ X 3 (r) = l + ^p. (36) 

The estimates X 2 (r), X 3 (r) are displayed in Figs. 2,3 (resp. 4,5), for comparison 
with the measured values X 2 (r), X 3 (r) at T/T c = 1.97 (resp. 5.7), i.e. (3 3 = 29 
(resp. 84). We see that the corrections to factorization implied by two-particle 
propagation provide a reasonable explanation of the behaviour observed for the 
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X : s at intermediate and short distances. This is especially true in the case of X2, 
showing that there is very little room for contributions from direct non-quadratic 
couplings in A 2 (x) in the full effective action (that resulting from integration over 
the gauge fields). This justifies our statement that the residual interactions between 
the colourless boundstates of the adjoint scalars are very weak. 



5 Conclusions 

In this paper, we have studied properties of the two dimensional model derived 
in by dimensional reduction of 3D QCD at high temperature. In this model, 
scalar fields A in the adjoint representation of SU(3) interacts via SU(3) gauge fields 
U, in addition to a self-interaction generated by integration over the non-static 3D 
gauge degrees of freedom. Such properties are interesting since it was shown in 



IIJ that dimensional reduction works remarkably well in this case. Also, they offer 
an opportunity to explore non-perturbative features in a low dimensional situation 
where the IR singularities are particularly severe. 

By means of numerical simulations, we explored that part of phase space where 
the i? T -symmetry A — > —A is unbroken, in accordance with the small A expansion 
used to derive the model, known to be valid quite soon above the transition tem- 
perature of pure 3D QCD. We identified two boundstates 5* and P, respectively 
even and odd under R T , and thus coupled to monomials respectively of degree 2n 
and 2n + 1 in the A's. The S signal coincides with that previously obtained from 



Polyakov loops correlations JTTJ] , where however the P-state contributions could not 
be disentangled. 

These results came out from the measurement of three even-even and three odd- 
odd distinct correlations, as functions of the on-axis lattice distance r. Great care 
was taken in the analysis of their shape in r, with the result that in all cases, the 
signal found was that expected from the occurrence of genuine poles in momentum 
space. A contrario, this demonstrates that the picture where the decay with r of such 
correlations reflects the propagation in 3D of p = 2n or 2n + 1 "electric gluons" , i.e. 
a correlation length equal to 1 jp times the "Debye screening length" , is inadequate 
in the case under study. 

By comparing the size of the three different correlations measured for each of 
the S and P sectors, we were able to show that residue factorization holds, as 
expected on general grounds when one particle propagates between different states. 
The agreement with factorization was expectedly found to be particularly good at 
large distances, but we could even show that deviations at shorter distances are 
to a large extent compatible with propagation of two particles, namely two S or 
S and P respectively in the S or P channel. The overall picture thus is that the 
scalar sector of the reduced model at large distances, thought to accurately describe 
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static properties of 3D QCD at high temperature, consists of two weakly interacting 
colourless particles, respectively even and odd under the R T symmetry of the model. 

There are several problems, which this study invites to investigate further. Of 
course, similarly detailed analysis for full QCD in (3+l)D would be interesting. 
Furthermore, the construction of a reduced model where the ^-symmetry of the 



pure gauge theory is not spoiled by the reduction process is highly desirable [27 



with the hope that it exhibits a transition to a symmetric Z$ phase analogous to 
the low temperature QCD phase. 
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Appendix: Mean Field Technique for Composite Fields 

Correlators 



A Formulae for Traces and Determinants 



Let be a complex N x N matrix. Here we give an elegant trick [gfj to compute 
the traces 

(Pn=tl^ (37) 

for n > N, given (j) p for p < N. 

Consider the determinant P/v(t) = Det(l — t<f>), where t is a complex variable. 
It is a polynomial of degree N in t and its term of degree N is (— l) N Det(4>), and 
we have 

log(Pjv(t)) = trlog(l-ty). (38) 

Both sides of this identity can be expanded in t in some finite neighbourhood of 
zero. The method consists in identifying the coefficients of the two series. The N 
first orders determine the coefficients of P/v(i) from the n 's ,n < N. Then, the 
higher orders directly express any <f> n , n > N as a function of the p 's, p < N. 
Note that instead of computing Det{4>) from the order iV coefficient of P^, one 
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can alternatively compute 0jv given Det((fi), which is convenient for SU(N) group 
matrices. 



If applied to <fi = A, an element of the SU(3) algebra, (in which case Ai — 0), 
this technique gives A n , n > 3 in terms of A 2 and A 3 taken as independent variables. 
The first non trivial identities are 

A, = \{A 2 )\ 

A 5 = ^A 2 A 3 , (39) 



1, . x, 1 

4— ' 3 : ' ,:! 



A 6 = -(A 2 f + -{A^ 2 



and the determinant is 



A, = \(A 2 y + ^(A 3 fA 2 , 
9 A 

DetA = (40) 



In what follows, we will have to manipulate monomials of the elementary scalar 
fields A a (x), defined for SU(N) through 

N 2 -l 

A(x) = A a {x)X a (41) 

a=l 

where the traceless basis \ a is subject to the normalization 

tr\ a \P = U aP . (42) 
On this basis the anti-commutators read 

N 2 -l 

X 13 } = c aP l N + d ^ x ^ ( 43 ) 

7=1 

with real and totally symmetric tensors c and d. With these normalizations and 
notations, we have 

N 2 -l 1 N 2 -l 

HA 2 = ^tY[\ a \P]A a AP = 5 a pA a AP, (44) 

a/3=l aP=l 

N 2 -l N 2 -l 

tiA 3 = tr [A^A^A 7 ] A a A l3 A 1 = - ^ d a(3l A a A p A\ (45) 

a/?7=l «/37=l 

16 



The projection property 



Af2_1 1 / 1 \ 

E = o (Ww-liFWcd) (46) 

a=l ^ ' 

can be used to derive that for any pair X, Y of complex N x N matrices the 
following identities hold: 

^tr\ a Xtr\ a Y = \ UtXY - ^trXtrA (47) 
EtrXA a FA a = i ftrXtrF - ^trXF j . (48) 

OL 

In what follows we specialize to N = 3. 



B Correlators of Composite Operators 

By gauge invariance, A 2 (x) and A^{x) can be chosen as the two effective degrees 
of freedom. By _R T -symmetry, the even and odd sectors under A — > —A decouple. 
Here we derive consequences of the assumption that their dynamics is determined 
at leading order by their given vacuum expectation values A 2 and respectively and 
their connected two-body correlations ^2,2(2^) and A^^x). 

We use a Wick like treatment to express the higher order connected correlation 
functions and averages through the quantities mentioned above. If n = 2p, any 
A a (x) is assigned to belong to a pair A 2 (x), then considered as a free field denoted 
S(x). So each monomial is replaced by a sum over all such pairings, and each of the 
p pairs is subject to the substitution W 2 , 

W 2 : A a A^5 a ^S(x), (49) 

leading to a monomial of degree p in S(x). If n — 2p + 1, one first performs the 
p — 1 possible substitutions W 2 (the result after p substitutions would transform 
as an octet under SU(3) and thus vanishes), which yield a monomial necessarily 
proportional to S p ~ 1 (x) tr A 3 (x). There we apply the substitution W3, 

W 3 : trA 3 (x) -> P(x), (50) 
where P(x) is also considered as a free field. 
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Once the local operators have been expressed in terms of S(x) and P(x), any 
average is obtained by using 



(S(x)) 
(S(x)S(0)} 
(P(x)S(O)} 

(P(x)) 



A 2 , 

A 2,2( X ) 

o, 
o. 



A 2 



(P(x)P(O)) = A 3i3 (x). 



(51) 
(52) 
(53) 
(54) 
(55) 



These rules generalize the way how in Section 4 we computed averages involving 
A^(x). Let us now detail calculations involving A 5 (x). 



From Eqs.(|M|, |4], |45|), we have 



6 



(56) 



We apply rule W 2 to the right hand side. The contraction of a with (3 produces 
S(x) P(x) once. Using Eqs. ([I7|, [|^, ^), one finds that the 6 W 2 — contractions 
of either a or (3 with either one of the three other indices contribute each the same 
amount 

\s{x) x UiA^x), 

that is from rule W3 

-S(x)P(x). {51 



(57) 



We thus arrive to the substitution 

A 5 (x) 



5 fi ^5 
^(l + *-)S(x)P(x) = ^S(x)P(x), 



(59) 



which we perform in the two point correlations A 3i5 (x) = (^3(0;) A 5 (0)) and ^5(0;) = 
(A 5 (x) A B {0)) to get 



A 5 , 5 {x) 



35 

-(P(x)P(0)S(0)), 



35y 
247 



(P(x)P(0)S(x)S(0)). 



According to Eqs. ( pl| - |55|) , these averages are given by 



35 



) A 3 , 3 (x) [A 2 2 + A 2)2 (x) 
18 



(60) 
(61) 

(62) 
(63) 



As a last application, we derive the value of ^3(0) = (A^(x)). By definition 



A 2 3 (x) = A a A (i A~ / ti\ a X l3 X y ^ 5 ^ ? trA 5 A £ A c , 

a /3 7 S e £ 



(64) 



where all the fields are taken at the same point x. Applying all possible W2 substi- 

(65) 



tutions and using Eqs. ([47], H) leads to the substitution 

5S 3 {x) 



A 2 Jx) 



64 



and averaging via Eq. (|5T| ) provides the final result ( |29| 
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Figure 1: The on-axis correlations A nm (r) at T/T c = 1.97 (f3 3 = 29), versus the 
distance in units of \/T. The even cases [n,m] = [2,2], [2,4] and [4,4] all have the 
same shape, and the odd cases [3,3], [3,5] and [5,5], again similar with each other 
in shape, are steeper. 
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Figure 2: Residue factorization: data for the quantity X2, Eq. (PH) at T/T c = 1.97 
(/?3 = 29) versus the distance in units of 1/T. It approaches one at large distances. 
The quantity X 2 corresponds to our interpretation (Section 4, Eq. (p5|)) of the 
deviation from one of X2 at shorter distances. 
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Figure 3: Residue factorization: data for the quantity X 3 , Eq. (|T^) T/T c = 1.97 
(/3 3 = 29) versus the distance in units of 1/T. It approaches one at large distances. 
The quantity X 3 corresponds to our interpretation (Section 4, Eq. (j36|)) of the 
deviation from one of X 3 at shorter distances. 
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Figure 4: Same as in Fig. 2 at T/T c = 5.7 = 84). 
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Figure 5: Same as in Fig. 3 at T/T c = 5.7 = 84). 
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Figure 6: The residues Ri defined by Eq. fl30|) stay close to one in the whole 
temperature range. 
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Figure 7: Plot of 4A24/5A2A22 at T/T c = 1.97 = 29). This quantity is one if 
Eq. (|3lD exactly holds. 
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